Anchoring of polymers by traps randomly placed on a line 
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We study dynamics of a Rouse polymer chain, which diffuses in a three- 
dimensional space under the constraint that one of its ends, referred to as the 
slip-link, may move only along a one-dimensional line containing randomly 
placed, immobile, perfect traps. For such a model we compute exactly the 
time evolution of the probability P s \(t) that the chain slip-link will not en- 
counter any of the traps until time t and consequently, that until this time 
the chain will remain mobile. 
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I. INTRODUCTION 



The time evolution of the survival probability P(t) of particles diffusing in a d- 
dimensional space in the presence of immobile, randomly placed traps has been widely 
discussed in the physical and mathematical literature within the last two decades. An in- 
terest in this problem has been inspired by the evident physical significance of the subject 
(excitation and charge motion, photoconductivity, photosynthesis). Further on, such an 
interest has been stimulated by an important observation [ffj that P(t) exhibits a non- mean- 
field long-time behavior, which is intimately related to the so-called Lifschitz singularities 
near the edge of the band in the density of states of a particle in quantum Lorentz gas 
and is reflected in the moment generating function of the so-called Wiener sausage 0,f|. 
Later works (see, e.g., have also pointed out the relevance of the issue to percolation, 

self-avoiding random walks or self-attracting polymers, as well as to the anomalous behav- 
ior of the ground-state energy of a Witten's toy Hamiltonian in supersymmetric quantum 
mechanics 0. 

Various analytical techniques have been elaborated to calculate P(t), including an ex- 
tension of the "optimal fluctuation" method pj, different methods of evaluating upper and 
lower bounds (see, e.g., @J3>0~il)> Green functions approaches ||, field-theoretic treatments 
||, as well as a variety of mean-field-type descriptions (see []T~D| |T^j| and references therein). 
These studies have revealed a two-stage decay pattern of the form 

, p ,.x f -rhrMt), t m <t<tc, (A) 

^ " \-n^h^\t>>t Cj (B) (1) 

where n tr denotes the mean density of traps, t m is a microscopic time scale and t c denotes the 
crossover time between the intermediate- (A) and long-time (B) kinetic stages. Further on, 
the function 4>d{t) appearing in Eq.(^.A) defines the mean volume of the so-called Wiener 
sausage (see, e.g. 0) — i.e. the mean volume swept by a diffusive spherical particle during 
time t. Its discrete-space counterpart, i.e. an analog of <pd{t) defined for lattice random 
walks, is referred to as the mean number of distinct sites visited by a particle up to the time 
t (see |16[] for more details). The functional form of (j>d{t) is different for different spatial 



dimensions d and obeys: 

t 1 ' 2 , d = 1 

o,/(f) :v [ t/ln(t), d = 2 (2) 
t, d>3 

The physical behavior underlying the kinetic regimes described by Eqs.(|l].A) and (]1|.B) 
has also been elucidated. It has been understood that Eqs.(jl].A) and (|l].B) are supported by 
completely different realizations of random walk trajectories: The intermediate-time behav- 
ior described by Eq.(^.A) is associated with typical realizations of random walk trajectories 
and is consistent with the predictions of the mean-field, Smoluchowski-type approaches 
|jH^|l§{j]|; namely, <j> d {t) = J* dt'k 

Smoi(^) @, where fcsmoi(^) is the so-called Smolu- 



chowski constant, which equals the diffusive current through the surface of an immobile 
(i-dimensional sphere. On the other hand, the long-time asymptotical form in Eq.([l|.B) 
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showing a slower time-dependence compared to the intermediate-time decay law, stems 
from the interplay between fluctuations in the spatial distribution of traps (namely, on the 
existence of rare but sufficiently large trap-free cavities), and atypical realizations of random 
walk trajectories which do not leave such cavities during the time of observation. Note also 
that Eq.(|I[B) describes the anomalous long-time tail of the moment generating function for 
the Wiener sausage volume @. 

Due to the general interest in fractal structures as useful approximate models of disor- 
dered media, and/or anomalous diffusion, Eq.([l[) has been extended |T^JT^ -[T6| to describe 



trapping of random walkers on random or regular structures characterized, in the general 
case, by a non- integer spatial dimension df and anomalous diffusion exponent d w , the latter 
being defined through the relation describing the time-dependence of the second moment of 
the particle's displacement, r 2 {t) ~ t 2 ^ dw , where d u may be different from the value d u = 2, 



which holds for conventional diffusive motion in Euclidean <i-dimensional space |14Hl6|. For 
such systems heuristic agruments []T2"| -fH|| suggest that P(t) follows an asymptotic behavior 
of the form 

, ■ -ThrMt), t m <t<t c , (A) 

lll^{t) OC < dL/Uf+cL) .A.KA.^A \ . . / 



with 



t d f / du ,^ dfKdu 
4>d(t) oc { t/ln(t), d f = d UJ (4) 
t, df > d^ 

Note, that here df can attain integer values and d w may be set equal to 2, which leads to 
conventional diffusive motion in Euclidean space; then Eq.(|[) reduces to Eq. ([]]). 

The decay patterns as in Eqs. (§) and @) have been verified numerically for different types 
of fractal systems, such as, e.g., Sierpinski gaskets or percolation clusters p|-|T6|j. However, 



rigorous results describing the evolution of P(t) in systems showing anomalous diffusion are 
lacking at present. 

In this paper we discuss a particular case of trapping; we consider reactions in systems 
in which the random motion of a particle moving in presence of randomly placed, immo- 
bile traps is non-Markovian, but, nonetheless, its survival probability can be determined 
exactly starting from first principles. More specifically, we consider a situation in which a 
mobile particle (referred to in what follows as a slip-link) is located at the end of a long 
polymer chain diffusing in three-dimensional space and is constrained to move along a one- 
dimensional line containing randomly placed perfect traps ([I]). The chain is modeled as a 
sequence of iV beads connected by phantom elastic springs and its dynamics is described 
within the framework of the customary Rouse model |T3 . 
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FIG. 1. Polymer chain with an active particle attached to one of its extremities (the slip-link). 



The paper is structured as follows: In Section 2 we formulate the model. In Section 3 
we discuss briefly the exact solution for the survival of a particle which diffuses on a line 
in the presence of randomly distributed traps, and then rederive this solution in terms of a 
path-integral method. In Section 4 we describe, following [19| the path integral formalism 
for evaluating the measure of trajectories covered by a tagged bead of a Rouse chain. Next, 
we show how such a formalism can be applied for the exact computation of the probability 
that the slip-link is not trapped until time t. Finally, in Section 5 we conclude with a 
summary and discussion of the obtained results. 



II. THE MODEL 

Consider a polymer chain embedded in three-dimensional space and consisting of iV 
beads (Fig.|TJ) , which are connected sequentially by phantom harmonic springs of rigidity K. 
The rigidity can be also expressed as K = 3T/6 2 , where T is the temperature (written in 
units of the Boltzmann constant ks) and b is the mean equilibrium distance between beads. 
Further more, the positions of all beads are denoted by fj = (xj,yj, Zj), where the subscript 
j enumerates the beads along the chain, j = [0, N]. All beads, except the slip-link (j = 0), 
may move freely in 3d. On contrary, we stipulate that the slip-link is constrained to move 
only along the X-axis, such that its position in space is given solely by the X-component 
of the vector r , since y and z must always equal zero. Then the potential energy t/({r}}) 
of such a chain is given by 



K N_1 



K 

1 - Ti ) = — 



j=0 



'JV-1 2 JV-1 

E ( X 3+X- X i) + E {(%'+! 
j=0 j=l 



Vj) 2 + ( z i+i - z jf} 
(5) 
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Next, apart from the holonomic constraints imposed by the springs, the beads experience 
the action of random forces fj(t). These random forces are assumed to be Gaussian and 
uncorrelated in time and space; their Cartesian components fj >a (t), where a = x,y,z, obey: 



fj,a{t)fj',a'{t') = 1C,T5jji5 a>a >5{t - t'), 

where the bar stands for averaging over thermal histories and ( is the macroscopic friction 
coefficient. 

In the absence of excluded-volume effects, the dynamics of the Rouse chain is guided by 
the corresponding Langevin-Rouse equations [|T8| 



sulfa}] 

C?i = — ^ + ( 7 ) 

where the dot denotes the time derivative. As one can verify readily, for the potential energy 
given by Eq. (||) the equations for the fj decouple with respect to the Cartesian components. 
That is, the dynamics of, independent of yj and Zj and reads 

(Xj = K(x j+1 + Xj-i - 2xj) + fj, x {t), (8) 

which equation holds for j = 1,...,N — 1. On the other hand, the displacements of the 
chain's extremities along the X-axis, i.e. x and xn, obey 

(x = Kfa - x Q ) + f , x (t), (9) 

and 

(x N = K(x N ^ - x N ) + f N>x (t) (10) 

Hence, with regards to polymer dynamics, one faces an effectively one-dimensional model. 

Lastly, we suppose that the iV-axis contains perfect, immobile traps, which are placed 
at random positions with mean density n tr . The positions of the traps are denoted by {X n }, 
—oo < n < oo. According to our model, the traps influence the dynamics of the polymer 
only by trapping (immobilising) the slip-link at the encounter. The influence of the trap 
on the other beads ( i.e. such that j e [1>-W]) is indirect: as soon as the slip-link gets 
immobilised by any of the traps, the chain becomes anchored as a whole due to the links 
between the chain's monomers. Our aim is to compute exactly the time evolution of the 
probability, P s \{t), that a polymer chain with a slip-link sliding along the X-axis will not 
encounter any of the traps (and thus will remain mobile until time t). 



III. MONOMER TRAPPING ON A LINE 

It seems instructive to recall first the time evolution of the survival probability in the 
simplest case, for X = 0, i.e. for a single chemically active monomer diffusing on a one- 
dimensional line and reacting with randomly placed, immobile traps. It is intuitively clear 
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that for a Rouse chain containing N beads, we should recover (apart from the renormalization 
of the diffusion coefficient) at sufficiently long times the behavior predicted for a single 
monomer, because for a finite chain the random motion of any of the chain's beads ultimately 



follows the conventional diffusion of the chain's center-of-mass ||18|| . On the other hand, 
at shorter times substantial deviations between the motion of the monomer and of the 
center-of-mass are to be found, due to the essentially non-diffusive characted of the slip-link 
motion, induced by the internal degrees of freedom of the polymer. This anomalous regime 
stemming from the internal relaxation modes of the chain will cause significant departures 
from the decay forms described by Eq. QT|) . The derivation of decay laws associated with this 
anomalous regime is the primary goal of the present paper and will be discussed in the next 
Section. 

As one may expect, in one-dimensional systems the situation simplifies considerably, 
since here the diffusive particle can react only with two neighboring traps and thus can- 
not leave the intertrap interval. This renders the problem exactly solvable, reducing it to 
the analysis of the particle survival inside a finite interval, followed by averaging over the 
distribution of the intertrap intervals. We outline such a calculation following the seminal 
method of Ref. M. 



A. The one-dimensional exact solution. 

We calculate first the probability *f>(x,x(0),t\L) that a diffusive particle (whose diffusion 
coefficient is D = T/Q, which starts at x(0) will not encounter the traps at Xo = and at 
X\ = L until time t. This probability follows as the solution of the following one-dimensional 
boundary problem: 

f d$(x,x(Q),t\L) _ D d 2 ty(x,x(0),t\L) 

Ot Ox 2 
*(x,x(0),t\L)\ x=0 = *(x,x(0),t\L)\ x = L = ( n ) 

k ${x, x(Q), t = 0\L) = 6{x - z(0)) 

The solution of Eqs.(|ll|) can be readily found by standard means and takes the form of a 
Fourier series: 

, v I v 2 ~ / 7i 2 n 2 Dt\ . /nnx\ . fimx(0)\ 
*(z,x(0),t|L) = - £ exp f —j sin (— j sin f — ±± \ , (12) 

Now, to compute the monomer survival probability we turn to the position-averaged function 

ty(t\L) = - / / dx(0)dxty(x,x(0),t\L), (13) 
L Jo Jo 

which can be computed from Eq.(|T2|) and reads 

mL) ^ g (S + !)-. exp {_^±^ (14) 

Next, the desired survival probability of the diffusive monomer P(t) is determined as the 
convolution 
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P(t) = / dL^ mon (t\L)V(L) (15) 
Jo 

where V(L) is the probability density of having a trap-free void of length L. For a completely 
random (Poisson) placement of traps V(L) reads: 

ViL) = n tT exp(-n tr L), (16) 

Consequently, one finds the following general expression determining the monomer survival 
probability 

P(t) = £ £ (2i + 1} - jT iL exp (- ^^ - ntli ) (17) 



The asymptotical behavior of the expression in Eq . (|17D has been discussed in detail in Ref. 
|J; it has been shown that P(t) follows a two-stage decay pattern as in Eq.(l). Explicitly 
on has 



/ exp (-AnUDt/Tr) 1 / 2 ) , t m < t < t c , (A) 

■ ~ » exp(-3(7rX^/4) 1/3 ), t > t c , (5) 



where the crossover time t c separating two regimes obeys t c ~ l/Dn 2 T , and consequently, 
can be large if n tr is small. 



B. Path integral solution of the monomer trapping problem in ID. 

In this subsection we rederive the solution of the monomer trapping problem in terms of 
the path-integral formalism, which will be later used to determine the trapping kinetics of 
the slip-link. To do this, we will proceed as follows: 

We first write the solution of Eq.flllf) as an integral over Brownian paths x(r): 

where V{x(t)} denotes integration over the monomer trajectories x(r), the exponential is 
the standard Wiener measure, while the subscript < x(r) < L, r G [0, t] signifies that the 
integral has to be calculated under the constraint that neither of the monomer's trajectories 
leaves the interval [0,L] within the time period [0,t]. 

Next, such a constraint can be automatically taken into account if we multiply the 
integrand by the step function depicted in the figure below 





1 


, U{X) 













(2m-l)L 2mL (2m + 1)L x 



(19) 
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Using the contour integral representation of such a step function [p0| 



1 r d\ e- Xx 



U{X) = 2ri I W 



1 if 2mL < x < (2m + l)L 



A(l + e Ai ) [ o if (2m + 1)L < x < (2m + 2)L 
where m = 0, 1, 2, . . ., we can rewrite Eq.(|T5D as 

1 



(20) 



$?(x,x(0),t\L 



, y 



- 1 C T dX e^(°) 



m=—oo 



A(l + e 



AM 



D{x(r)}exp -S{x(r)} , (21) 



where the action S{x(t)} is given by 



s(0) 



+ A 



dx(r) 
dr 



(22) 



Now, to compute the path-integral in Eq.([H]) with the quadratic action in Eq . (p22|) 
we have merely to define the action-minimizing trajectory x(r) and calculate the action 
corresponding to such a trajectory. The action-minimizing trajectory is defined by the 
classical Euler equation of motion, which for the action in Eq.([22|) is simply 



d f 1 dx(r) 



dr \ 2D dr 



A = 



Integrating this equation subject to the conditions x(r = 0) = x(0) and x{j = t) = x, we 
find 



T 



x( T ) = (x-x(0))- + x(0) 



(23) 



Consequently, the minimal action is given by 



and hence, the formal solution of the boundary problem in Eq.([TT]) can be written down as 

1 ^ 1 c+ l°° j\ r r^_wn^ 2 i 

tf(z,z(0),t|L) 

(24) 



J oo 

2 ex p 



2vrz 7 A(l + e XL ) " 1 4Dt 

c— ioo 

r (x -a;(0) -2mL) 2 l 



Next, using the well-known representation of Jacobi theta-function 2^ 

g q (^) 2 = 9 3 La, exp (^W In" 1 / 2 I 

and setting 
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S 



exp 




x — x(0) 
2L 



(25) 



we may rewrite Eq.(f2|) as 



2 ~ / Vn 2 m 2 t\ f irm(x - x(0)Y 
*(x,x(0),t\L) = - J2 ex P JjT- cos L 



m=—oo 



( 




m=—oo 



Dn 2 m 2 t 
T 2 



/ 7imx\ { irmx(0)\ ( nmx\ ( 7rmx(0) 

sin I — ; — ) sin | ; + cos ( — ; — ) cos 



V L 



V 



V L 



\ L 



\ L 



(26) 



part A 



part B 



/ 



Note now that the contribution "A" (a product of two "odd" sine-functions — see Eg. (|26|) ) 
vanishes at x = and x = L. It thus mirrors the absorbing boundary conditions rele- 
vant to the trapping problem, while the second contribution "B" (a product of tho "even" 
cosine-functions) corresponds to the totally reflecting boundary conditions. Moreover, the 
expression in Eq.(^) in which only the odd contribution "A" is taken into account coin- 
cides with the result presented in Eq.([L2]). It follows that the asymptotical behavior of the 
monomer survival probability can be evaluated in terms of the path-integral formalism by 
picking then the odd contribution. We would like to note also that the construction we have 
employed is nothing but the well known "mirror principle", widely used in electrostatics. 

We conclude this subsection with the following prescription: 
Prescription (the "mirror principle"). To find the solution of the diffusion problem 
with the Dirichlet boundary condition at the ends of the segment [0, L], we have merely to: 
a) obtain the Green's function solution ^(x, x(0), t) on the full line {x,x(0)} g] — oo,oo[; b) 
restrict x — x(0) to the segment [0, L], perform the replacement x — x(0) —>■ x — x(0) — 2mL 
and take the odd contribution to the sum J2m=-oo ~ x (0) — 2mL,t). 



IV. TRAPPING OF THE SLIP LINK OF A POLYMER 



To calculate the evolution of P s i(^) in the case of a slip-link attached to a polymer 
chain, we will proceed essentially along the lines of the previous section. First, we present 
the derivation of \I/ s i(xo, x (0), tj — the probability distribution for the displacements of the 
slip-link on an infinite line without traps . Then, using the above formulated prescription 
(the "mirror principle"), we determine ^si^Xo, ^o(O), t\Lj — the probability that the slip-link 
of the chain, which is initially located at some point xq(0) inside the interval [0, L\, will not 
leave this interval until time t. Finally, the desired probability P s \(t) will be obtained from 
^/ s i(x , ^o(O), t\L) by averaging over the Poisson distribution of the interval's lengths. 



A. The probability distribution of the slip link displacement 



In this subsection we outline, following the analysis of Ref. ||19|| , the steps involved in 
the derivation of the probability distribution ^si^xo, Xq(0), tj for the displacements of the 
slip-link on an infinite line. 
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This probability can be written as 



*si(a;o,Xo(0),t) = f ... f dx 1 ...dx N ^(x ,x 1 ..,x N ,t) (27) 

v ' J—oo J—oo 

where {xo, xi, ■■, %n} = {xi(t), ..,XN(t)} here and henceforth denote the coordinates of all 
chain segments at time t, while ^(xo,xi, ..,xjsr,t) is the joint distribution function of the N 
beads of the polymer chain. In other words, ^/(xq,Xi, ...,x^,t) is the probability of having 
the x-coordinates of the N beads of the chain at the positions {xj(t)}, provided that initially 
they where at {xj(0)}. 

For the Rouse chain whose potential energy obeys Eq.fl^), the time evolution of the func- 
tion ty(xo, Xx, .., Xn, t) is governed by the following Smoluchowski-Fokker-Planck equation 



d N d ( d d 

-*(*o, x 1; .., x N , t) = D g — y— + P^-U(x , 



xi,..,x N ) )^(x ,x 1 ,..,x N ,t) (28) 



where (3 = 1/T and D = T/( is the diffusion constant of an individual bead (monomer). 
Equation (|28D has to be solved subject to the initial condition: 

N 

*{x o ,x 1 ,..,x N ,t = 0) = H5(x j (t)-x j (0)) (29) 

3=0 

Eqs . (|28|) — (|29|) define completely the evolution function ^(x , xi, ..,x^,t). 

The computation of the probability distribution ty(x o ,x o (0),t) using the path-integral 
formalism was first performed in [pl|] . Let us recall the main steps of this approach. First of 
all, it is expedient to cast Eq . (|28|) into the form of the (A^+ 1) -dimensional Schrodinger-type 
equation. This can be readily performed by making use the following ansatz: 

^(x ,xi, ..,x N ,t) = exp — t ~ 7} U ( x o, x h ..,Xiv)| $(x ,xi, ..,x N ,t), (30) 

where the function $(x ,^i, ..,XN,t) satisfies 

d N ( d 2 \ 

— $(x , xi, .., x N , t) = ~ P 2k2 ( x j-i ~ 2x i + x i+i) 2 J ^(^o, x N , t) (31) 



Now, by virtue of the Feynmann-Kac theorem, the formal solution of Eq. (|3l| ) can be written 
down explicitly as the following path-integral 



1 XO XN N 

$(x ,x 1 ,..,x N ,t) = jj : J ... J J[ P{a; i (r)}exp [ - S {x (r) , Xx(t) , ..,x n (t)}] (32) 

xo(0) Xiv(0) j=0 

where J\f is the normalisation constant and the action S^Xq(t), Xi(t), .., Xiv(r)| has the 
form: 
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S{x (t),xi(t), ..,x n (t)} 



dr 



1_N(dXj(T) 

AD ^ 

^ u 3=0 



N-1 



dr 



Df3 2 K 2 E (^i-i(r) - 2 Xj (r) + x ]+1 {r)f 



(33) 



In the subsequent calculations we will not specify the normalisation constant and pre- 
exponential factors; all our results will hold thus only to exponential accuracy. 

Now, we suppose following Ref. ]19| that the slip-link moves along some prescribed 
trajectory X A { matter of fact, such a constraint allows to symmetrize the boundary 

conditions for the action-minimizing trajectory, which are otherwise different at different 
chain's extremities. Such a constraint can be taken into account by multiplying the integrand 
in Eq. (^) by a functional delta-function of the form (see Ref. |19| for details) 



6(x (t) - X a (t)) = I V{m{r)}expi-t f dr m(r)(x (r) - X A {r)) 



which means that the action in Eq.(^) is replaced by an effective action of the following 
form: 



S'{x (t),Xi(t), ...,x n (t)} 



1 



N 



drl—Y 
1 AD ^ 



3=0 



dxjij) 



+ 



N-1 



(34) 



D(3 2 K 2 E (x^) - 2x 3 (r) + x J+1 (r)) 2 + im(r)(x (r) - X A (rj) 



3=1 



and one has to perform afterwards an additional integration over the measure T>{m{r)}. 
Now, the action in ([32]) is minimal for classical trajectories satisfying the Euler equation: 

£ (xo,X!, ...,x N ,x ,xi, ..,x N ) = 0; j G [0, N] , 

where £ is the Lagrangian function: 

2 N-1 







- d ) 




{dxj) 


dxj) 



£ 



1 "(dxfr) 



AD 



3=0 



dr 



D(3 2 K 2 E - 2x j (r) + x j+1 (rj) + im(r)((x (r) - X sl (r] 

3=1 

(35) 



Turning to the continuous j-limit, one obtains the following Euler equation, determining 
aje 





the optimal trajectories of the chain's monomers |19 

Q2\ 



dr 



-wo- 



d d 2 \ 

+ 2K(3D—\x j {r)=ADi5{])m(T) 



8t ' dj 2 

which has to be solved subject to the boundary conditions (see |l 
dxj(r 



(36) 



dj 
dr 



0. 



d 3 Xj(r) 



j=0,N 



-2Kf3D 



dp 
8 2 xAt 







3=0,N 



df 



T=0 



dr 



df 



(37) 



T = t 



11 



The action-minimizing trajectories Xj(r), defined by the boundary problem (|36] 



are obtained explicitly |19| in form of a series expansion over the normal Rouse modes ||18|| : 



2iN ~ _ 2 / Trpj \ ft j \r-r\ 2 I 

: 2jp cos ) dr m(T ) exp j J ^ ' 



where tr = (N 2 /2n 2 K is the largest fundamental relaxation time of the harmonic chain, 
i.e., the so-called Rouse time. This time may be interpreted as being the time needed for 
some local defect, e.g., kink, to spread out diffusively along the arclength of the chain. 

Next, substituting the expression for the optimal trajectory in Eq.(^) into Eq.(|35|) and 
performing the integration over T>{m(r)}, one arrives at the following general result 

*..(-„ ( o), ( )4 7 ^>M-/:(^) *jf (^) « 



(r - t') 



(39) 



where <f>(r — r') is given by 



r — r 



1^-4 k-T'|>TH, (A) 

^ l < 2 W-r T V\\T-T>\<rn, (B) <40) 



1? 



Equations (|39|) and ([5]) represent the desired generalization of the classical Wiener result 
for the measure of Brownian particle trajectories to the more complicated case of a particle 
attached to a diffusive Rouse chain. 

Finally, in order to compute the probability \l/ s i(xo, xo(0), that the slip-link of the 
Rouse chain will remain until time t within the interval [0, L\, we make use of the "mirror 
principle" prescription of the previous section. Multiplying the integrand in Eq . (|39|) by a 
step-function, we get 

*^o,x o (0),t|L)=^ £ — / — T ^ I - / D{X sl (r)}exp[-4x sl (r)}] (41) 

m ~ 00 c-ioo x 

where 

- 1: (^) (^) ^ - «o + v> (42) 

Below we discuss the asymptotical forms of ty s i(xo,xo(0),t\L\ using Eqs. (f4~l"|) and (f42"D 
in the limits |r — r'\ < t r and |r — r'\ > t r . 



B. Asymptotic behavior of ^ s \[xq, xq(0), t\L J . 

We focus first on the behavior in the intermediate-time limit, t < tr. Note, however, 
that since tr ~ N 2 , for sufficiently long chains this intermediate-time regime may last over 
quite an extended time interval. 
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It follows from Eqs. (|41| ) and (|42|) that for t < tr the internal relaxations of the chain are 
most important and lead to the following form of the action 

- (£ Ti: m *r m ^ w> A 



(43) 

Now S in Eq.|43| is non-local and possesses non- Wiener scaling properties. In particular, it 
yields for the mean-square displacement of any chain's bead, including the slip-link, the 



following law fli3,18|,19 



^(t)-t 1 / 2 , (44) 

i.e. a subdiffusive behavior, which signifies that (for the time scales considered here) the 
trajectories of the chain's beads are spatially more confined that there of simple Brownian 
particle. 

Now, one can readily find that the optimal trajectory X s \(t) which minimizes the action 
in Eq.([J3|) obeys the following Euler equation 



We seek the solution of Eq.(^5|) in the form 

l sl (r') = A (rT 

where r' = ru and < u < 1. Substituting this form into Eq. (|4"5|) , we get the functional 
equation 

A^( l **L+X = C (46) 



\D ) Jo jT^u 

where C is some constant. Note now, that the left-hand side of Eq . (f46|) is independent of r 
for a = —1/2 only, which thus fixes the value of a to a = —1/2. Hence, we obtain 

x ^) = — bid (47) 



which leads to the following expression for the minimal action 

S{X sX {r)} = - U±j (C - A) 2 + A (x Q - * (0)) (48) 
In turn, the value of the constant C can be found by integrating Eq. (f47|) , which gives: 

1 /2 

C = X + l(jti) (*o-* (0)) 
13 



Substituting this expression into Eq. (fi8|) , we arrive at the final equation for the minimal 
action 

S{X sl (r)} = \ ' (xq - x (0)) 2 + A (x - xo(0)) (49) 

Next, taking advantage of Eqs. (|49|) and (|4lD , we find that the probability that the slip- 
link which is at point x at t = will stay inside the interval [0, L] until time t is given to 
exponential accuracy by 

# s i(zo,Xo(0),*|L) « £ ex p|"^(^) (x -x (0)-2mL) 2 |, (50) 
which yields, by virtue of the "mirror principle", the following result 

Further on, integrating Eq. (B3f) over x and a; (0), we get for the position-averaged function 
V A (t\L) (see Eq.(13)) 

and consequently, the desired probability P s \(t) that the slip-link will not encounter any of 
the traps until time t is given by the following integral 

00 r°° ( 7r 2 (21 + l) 2 \ 
P sl (t) « £(2Z + l)- 2 ^ rfL exp ^ -^-DO " ntrL j (53) 

where for notational convenience we have introduced the "effective" time 9, where 6 = 
(t/pKD) 1 / 2 /2n. Note that Eq.(g) becomes identical to Eq.fllTD upon the mere replacement 
£ — > 0, which readily enables us to get the corresponding decay forms from Eqs. (|T8D . We 
recall, however, that Eq.(^) is valid only for t < tr, which will result in a slightly more 
complicated overall decay pattern than the one described by Eqs. ([18]). 

Consider next the evolution of ^f s i(x , xq(0), t\Lj and, respectively, of P s \(t) in the limit 
t > tr. Note that here the action in Eq . (|42|) reduces to the standard result 



which is simply the action of an isolated Brownian particle which moves with the diffusion 
coefficient D/N (compare to Eq.(|2"2])). Consequently, in this time limit we have that the 
mean-square displacement of the slip-link obeys 

~xW)~^t (55) 
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and the probability P s \(t) follows 



, v 2 f 00 , ( 7T 2 (2/ + l) 2 Dt X 

P sl (t) « E( 2/ + !) y cLLexpj^ n tr L 



(56) 



Note that this result could be expected on intuitive grounds, since, as we have already 
mentioned, in the limit t > tr the motion of every bead of the chain follows mainly that of 
chain's center-of-mass. 



C. Trapping pattern P s \(t) for the slip link of a Rouse chain. 

Consider first the case of an infinitely long chain, then tr = oo and the dynamics of the 
slip-link is described by Eqs.(|39|),([40|.B) and (f44l) over the entire time domain. Comparing 
the decay forms in Eqs.(|I7]),(|]l|) and fl53|), we readily find that the probability that the 
slip-link will not encounter any of the traps until time t shows the following two stage decay 
pattern: 



exp (-^n tI UDt//3K) 1/ J , for t m < t < t Csl , (A) 

P(£) ~ | \ V 16 ffl) 

exp (-ln 2 J 3 (7r 2 Dt/(3K) 1/6 ), fort > t c>1 , (B) 



In Eq. (|57|) the crossover time t Cj i separating the regimes A and B obeys t Cj i ~ j3K/Dn^ T . 
Note that t Cj i ~ (bn tr )~ 2 t c , where b is the mean distance between the chain's beads in 
equilibrium and t c is the corresponding crossover time in Eq. (|i8|) . Hence, for systems with 
a small density of traps t Cj i can be significantly larger than the monomer crossover time 
t c in Eq. (|T8|) . Note also that the Eq. Q5"7|. A) is the mean-field, Smoluchowski-type result 
corresponding to the sub-diffusive motion of the slip-link, described by Eq.([H]), in presence 
of uniformly distributed traps; the exponent in Eq.Q5"7|.A) is just the product of the trap 
mean density and the mean maximal range (span) of the slip-link displacement. On the 
other hand, Eq. (|57| .B) stems from the interplay between exponentially rare, large trap-free 
voids, (the large-L tail of Eq.(J16|)), and anomalously confined trajectories of the slip-link. 
Equation (|57| .B) also describes the long-time tail of the moment generating function of the 
Wiener sausage volume for the the slip-link trajectories of an infinitely long chain. 

We turn next to the case of finite chains, and hence to finite r^, which sets the upper 
bound on the time of applicability of the decay pattern in Eq. (p?[) . Here, at times greater 
than tr the conventional diffusive motion of the slip-link is restored, Eqs.(|39|),(|40|.A) and 
fl55|), and the decay has a form similar to that in Eqs.(l), 

p( f exp (-AnUDt/nN) 1 / 2 ) , r R « t « t c , 2 , (A) 

S exp (s^nlDt/AN) 1 / 3 ) , t > f C)2 . (B) 



Here the crossover time t Cj 2 between the A and B regimes is given by t c ^ ~ NJ Dn 2 r = Nt c . 

Note, however, that the overall decay pattern of P s \(t) is not the does not necessarily 
follow sequentially after Eq. ([5TD , i.e. the sequence given by Eqs.fl5"T|.A), (|57|.B), (|58|.A) and 
finally, (|58| .B) may be realized only if the crossover times would obey the following multiple 
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inequality t m <C i c> \ <C tr <C t c ^ which practically is never the case. To show this explicitly 
and to construct the actual overall decay pattern, it is expedient to rewrite the crossover 
times t c ,i and t Cj2 in terms of the Rouse time r R . We have then £ Cjl « Q~ a t r and £ C)2 ~ Q~ 2 tr, 
where Q = n tr (WV 1//2 ), i.e. is equal to the mean number of traps in the area covered by a 
Rouse chain of arclength bN in its typical equilibrium configuration. 

Below we analyse different possible situations with respect to the values of the parameters 
Q and q = n tr b and discuss the corresponding decay patterns. 

Case I. High density of traps, q ~ 1, and long chains, Q « 1. Note first 
that in this case the crossover time t Cy \ is comparable to the microscopic diffusion time, i.e. 
t c ,i ~ t m ~ b 2 /D, which implies that the Smoluchowski-type regime corresponding to the 
sub-diffusive slip-link motion, Eq.(^7LA), is unobservable. Hence, the law in Eq.(|5T|.B) will 
describe the decay of P s \(t) until t ~ t r . Further on, the condition Q > 1 implies that 
t c ,2 <C tr and the regime described by Eq.([58|.A) does not exist. Consequently, the overall 
decay pattern in the case I reads 




3„ 



njf (ix 2 Dt/f3K) ) j t m t <C Tr, (A) 



/>,(/) « I ' V 2 v --•<--, j " (59) 

' -3n 2 f(n 2 Dt/4N)^) , t » r R , (B) 



Note that both A and B regimes are essentially non-mean-field and stem from the presence 
of fluctuation trap-free voids. We also remark that in Eqs.(|59|) the most representative 
regime is the one in Eq.(^.A), which is associated with the sub-diffusive behavior of the 
slip-link, Eq. (f44D . Note that the value of P s \(t) at the crossover time separating the A and 
B regimes, i.e. P s \(t = tr), is of order of exp(— Q 2 ^ 3 ) 1, which means that the probability 
that the slip-link will be trapped during the stage described by Eq.(^.A) is considerably 
higher than the probability that it will be trapped according to the law in Eq. (|59"|.B). 



Case II. Low density of traps, g« 1, and very long chains, Q > 1. In 
this case t Cj i 3> t m , which implies that the regime in Eq.(|57|.A) describes the initial kinetic 
stage. Next, since here t c i tr, the regime in Eq.(|57].B) will also exist and will describe 
the intermediate-time kinetic behavior of P s i(^)- Lastly, the final stage will follow the decay 
in Eq.(^.B), because t Cj2 appears to be much less than the Rouse time tr and thus the 
regime in Eq.(|58|.A) will be absent. Consequently, in case II one has the following overall 
decay pattern: 



exp (-^n tr {ADt/pK) 1/4 ^ , t m < t < t Cjl , (A) 
exp [-In^^Dt/PK) 1 ^ , t Cjl < t < r R , (B) (60) 
^ exp (s^nlDt/AN) 1 / 3 ) , t > r R , (C) 



In Eqs. ( |S0|) the decay laws in the first two lines are associated with the sub-diffusive motion 
of the slip-link and describe the mean-field, Smoluchowski-type (A) and fluctuation-induced 
(B) kinetic stages, respectively, while the law in the third line describes the survival of a 
Brownian particle (with diffusion coefficient D/N) in the fluctuation trap-free voids. Note 
that here the most representative regime is the one associated with the sub-diffusive motion 
in the slip-link and fluctuation trap-free voids, namely, regime (B). As in case I, we have 
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here P s \(t = tr) ~ exp(— Q 2 ^ 3 ), which is very small. This implies that the most probable 
decay is given by Eq.(^.B). The first stage, i.e. the decay described by Eq.(pD|.A), appears 
to be relatively unimportant, since during this stage P s i(t) does not drop appreciably, P s \(t = 
t c ,i) « 1/3. 

Case III. Low density of traps, g<l, and short chains, Q <C 1. In this par- 
ticular situation we have that t Cj i 3> t m and consequently, the initial decay obeys Eg. (|57|.A). 
Further on, since here t Cj i also exceeds the Rouse time, i.e. t Cj i 3> tr, the regime predicted 
by Eq.(|57|.B) is absent, which means that the decay in Eq.([57|.A) crosses over at t = tr to the 
decay predicted by Eq.(|58|.A). Lastly, the stretched-exponential dependence in Eq. (|58|.A) is 
followed at t > t c<2 , t c2 ^> tr by the form Eq.(^.B). Hence, in case III one has that P s \(t) 
follows 



exp(^-fn tr (4 J Dt/^) 1/4 ) , t m < t < t r , (A) 
exp(-4n tr (Dt/nNy/ 2 ), r R < t < t c>2 , (B) (61) 
exp (-Znll z {Tt 2 Dt/AN)y*) , t > t c>2 , (C) 



In this case, however, only the last, fluctuation-induced regime C appears to be significant; 
one can readily verify that P s i(^) practically does not change during the regimes A and B, 
P s l{t = L 2 ) w exp(-4/7r) ~ 1. 



V. CONCLUSIONS 

To summarize, we have studied the dynamics of an isolated Rouse chain, which diffuses 
in a three-dimensional space under the constraint that one of its extremities, the slip-link, 
may move only along a line containing randomly placed immobile traps. For such a model 
we have computed exactly the time evolution of the probability P s i(^) that the slip-link will 
not encounter any of the traps until time t, i.e. that the chain will remain completely mobile 
until this moment of time. We have shown that in the most general case this probability is a 
succession of several stretched-exponential functions of time, where the dynamical exponents 
depend on the time of observation and on characteristic crossover times. We have specified 
these crossover times and have determined explicitly the forms of P s \{t) in several particular 
situations. 
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